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SUMMARY 


Equations  of  motion  and  associated  boundary  conditions  are  developed  for 
the  general  nonlinear  vibrational  behavior  of  thin  conical  shells.  The 
theory  is  based  upon  nonlinear  strain-displacement  relations  deduced  for 
a  conical  shell  from  those  derived  by  Sanders  for  thin  shells  of  compound 
curvature.  Equations  for  the  bending,  buckling,  and  postbuckling  of 
conical  shells  under  arbitrary  loads  are  developed  also  and  are  shown  to 
reduce  to  equations  based  on  more  simplified  theories  for  both  conical 
and  circular  cylindrical  shells  and  circular  flat  plates.  Various 
solution  approaches  to  the  nonlinear  conical  shell  vibration  problem  are 
examined,  and  a  new  numerical  method  of  solution  is  proposed  and 
discussed. 


FOREWORD 


The  work  reported  herein  constitutes  a  portion  of  a  continuing  effort 
being  undertaken  at  Stanford  University  for  the  U.  S.  Army  Aviation 
Materiel  Laboratories  under  Contract  DAAJ02-68-C-0035  (Task  1F162204A17002) 
to  establish  accurate  theoretical  prediction  capability  for  the  static  and 
dynamic  behavior  of  aircraft  structural  components  utilizing  both  con¬ 
ventional  and  unconventional  materials.  Predecessor  contracts  supported 
investigations  which  led,  in  part,  to  the  results  presented  in  References 
18,  28,  29,  32,  and  33. 
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INTRODUCTION 


All  investigations  of  he  free  vibrations  of  conical  shells  to  date  have 
used  linear  theory,  mostly  based  upon  Love's  first  approximation  theory 

for  thin  shells*  .  Some  of  the  notable  earlier  studies  are  those  by 

2  14  5  6 

Strutt  ,  Federhofer  ,  Goldberg  ,  Grigolyuk  ,  Shuman  ,  Saunders,  Wisniewski, 

7  8 

and  Paslay  ,  and  Garnet  .  The  last-named  author  modified  the  theory  to 
take  into  account  transverse  shear  and  rotatory  inertia  effects. 

With  the  development  in  the  late  1950s  of  the  more  sophisticated  theory 

for  conical  shells  (developed  specifically  for  buckling  problems)  by 
9  10 

Seide  and  Singer  ,  free  vibration  studies  employing  this  theory  were 

11  12 

later  carried  out  by  Seide  and  Weingarten  .  These  studies  were  based 
on  the  linear  strain-displacement  relations  given  by  Love's  first  approxi¬ 
mation  with  circumferential  displacements  in  the  curvature  terms  neglected 
and  with  midsurface  inertia  terms  omitted.  Results  of  these  studies 
correlated  well  with  experimental  data  for  vibrating  conical  shells 

exhibiting  a  large  number  of  circumferential  waves.  This  was  expected, 

13 

since  neglect  of  the  midsurface  inertias  was  shown  by  Reissner  to  be 
valid  only  when  a  large  number  of  circumferential  waves  is  present. 

14  15 

In  two  other  studies,  notably  those  by  Chen  and  Weingarten  ,  the 

authors  retained  the  circumferential  terms  in  the  curvature  expressions 

and  also  included  the  previously  neglected  midsurface  inertia  terms. 

14 

Chen  used  a  Lagrange-equation  approach  in  conjunction  with  the 

Rayleigh-Ritz  method  for  a  conical  shell  with  classical  simple  supports. 

Again,  correlation  with  experimental  data  was  good  for  higher  mode 

numbers,  whereas  results  for  a  small  number  of  waves  correlated  poorly. 

15  18 

Weingarten  used  Saunders'  linear  shell  equations  in  conjunction  with 
a  finite-difference  scheme  introduced  by  Budiansky  and  Radkowski*^  for 
static  analysis  of  shells.  The  effects  of  different  boundary  conditions 
upon  the  vibrations  were  assessed,  and  it  was  shown  that  at  low  wave 
numbers  the  boundary  conditions  have  a  marked  effect  upon  the  frequencies. 
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This  was  also  shown  in  the  earlier  work  of  Seide^  for  two  different 
types  of  circumferential  boundary  restraints  of  simply  supported  conical 
che 11s. 

To  the  best  knowledge  of  the  authors,  no  investigation  of  the  vibrations 
of  thin  conical  shells  using  nonlinear  theory  has  been  carried  out  to 
date.  It  is  therefore  the  purpose  of  the  present  work  to  develop  the 
governing  equations  and  the  associated  boundary  conditions  for  the  non- 
liuear  vibrational  behavior  of  a  thin  truncated  conical  shell  and  to 
examine  and  discuss  various  possible  approaches  for  solving  these  more 
accurate  but  complex  equations. 

Recent  investigations  of  the  related  problem  of  nonlinear  vibrations  of  a 
circular  cylindrical  shell  have  produced  some  interesting  results.  In  a 

major  work  on  the  subject,  in  which  midsurface  inertia  terms  are  neglected, 

18  19  20  21  22 

Mayers  and  Wrenn  criticize  previous  solutions  *  *  ’  based  upon  the 

well-known  von  Karman-Donnell  strain-displacement  relations.  The  authors 

carry  out  a  new  solution  for  the  classical  equations  which  removes  the 

basis  for  criticism  and  discloses  the  existence  of  a  nonperiodic  vibration 

behavior,  a  phenomenon  not  considered  previously.  Solutions  for  a  small 

number  of  circumferential  waves  are  obtained  also  using  the  more  accurate 

23 

strain-displacement  relations  derived  by  Sanders 

In  the  present  invescigaticn,  the  governing  equations  for  the  nonlinear 

vibrational  behavior  of  a  conical  shell  are  developed  using  nonlinear 

strain-displacement  relations  for  a  conical  shell  deduced  from  those 
23 

derived  by  Sanders  for  a  shell  of  arbitrary  shape.  The  total  strain 
and  kinetic  energies  of  the  shell  are  developed  for  application  in 
Hamilton's  variational  principle;  the  governing  differential  equations  of 
motion  and  associated  boundary  conditions  are  obtained  therefrom.  As  a 
special  case,  these  equations  are  shown  to  reduce  to  the  so-called  "clas¬ 
sical"  equations  of  conical  shells;  namely,  those  for  which  the  circum¬ 
ferential  terms  in  the  curvature  and  nonlinear  midsurface  strain  expres¬ 
sions  are  omitted  and  the  midsurface  inertia  terms  are  neglected. 

Although  the  present  work  deals  with  the  nonlinear  vibrations  behavior 
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of  a  conical  shell,  with  little  added  effort  the  governing  equations 
and  associated  boundary  conditions  for  the  nonlinear  behavior  of  a 
conical  shell  under  arbitrary  static  and  dynamic  loads  are  obtained  also. 
These  equations  are  valid  for  either  forced  nonlinear  vibration  problems, 
nonlinear  vibration  of  a  prestressed  shell, or  prebuckli&g  and  postbuckling 
nonlinear  behavior  of  an  arbitrarily  loaded  shell.  These  extended 
equations  are  presented  in  Appendix  I.  Linearized  stability  equations 
and  boundary  conditions  using  the  nonlinear  terms  developed  in  this 
study  are  also  presented  in  Appendix  II. 

The  equations  derived  throughout  this  investigation  are  quite  general  in 
that  they  reduce  (by  a  set  of  simple  substitutions)  to  corresponding 
equations  for  circular  cylindrical  shells  and  circular  flat  plates.  These 
equations  are  presented  in  Appendixes  IH  and  IV,  respectively,  and  they  are 
compared  with  equations  derived  in  previous  studies. 

Lastly,  in  attempting  to  solve  the  nonlinear  vibration  problem  for  a 
conical  shell,  related  work  done  for  the  limiting  case  of  the  circular 
cylindrical  shell  is  examined. 

The  first  known  investigation  of  the  nonlinear  vibrations  of  a  circular 

24 

cylindrical  shell  was  carried  out  by  Reissner  using  the  classical 
equations  (i.e.,  with  midsurface  inertia  neglected).  For  linear 
vibrations,  Reissner  assumed  that  the  radial  displacement  function  may  be 
represented  as  a  chessboard  pattern.  This  assumption  was  validated  by 
previous  and  xecent  experimental  work  as  reviewed  in  Reference  18.  For 
nonlinear  vibrations,  Reissner  assumed  that  the  nonlinearity  has  a  more 
pronounced  effect  on  the  arbitrary  time  function,  which  modifies  the 
choice  of  deflected  shape,  than  on  the  deflected  shape  itself;  conse¬ 
quently,  Reissner  also  used  the  chessboard  pattern  for  the  nonlinear 

21 

vibration  solution.  However,  as  pointed  out  by  Evensen  ,  the  use  of  the 
chessboard  pattern  leads  to  a  circumferential  displacement  that  is  not  a 
periodic  function  of  the  circumferential  coordinate. 

Another  solution,  which  like  Reissner'?  did  not  satisfy  the  periodicity 

19  20 

condition,  was  presented  by  Chu  .  In  a  later  study  by  Nowinski  ,  the 
periodicity  was  taken  into  account;  but  as  a  consequence  (as  noted  by 
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Evensen  ),  a  nonzero  radial  displacement  occurred  at  the  end  of  the 

shell,  thus  violating  the  assumed  geometric  boundary  conditions  of  free 

support.  All  results  obtained  in  the  foregoing  references  vere  based 

18 

upon  periodic  motion,  except  those  of  Mayers  and  Wrenn  . 

The  nonlinear  behavior  experienced  by  a  circular  cylindrical  shell  should 

be  experienced  as  well  (perhaps  in  a  slightly  modified  form)  by  a  conical 

shell.  Consequently,  in  SUGGESTED  METHODS  OF  SOLUTION,  extension  to  a 

conical  shell  of  the  methods  of  solution  used  by  the  cylindrical  shell 

investigators  is  discussed  and  an  alternative  solution  based  upon  a 

modified-Reissner  variational  principle  approach  (as  introduced  by  Mayers 
28  29 

et  al.  ’  )  is  also  suggested.  Finally,  a  new  direct  numerical  scheme 

for  solution  is  proposed  and  demonstrated. 


STRAIN-DISPLACEMENT  REIATIONS 

The  nonlinear  strain-displacement  relations  employed  in  this  study  are 

23 

deduced  from  those  derived  by  Sanders  for  thin  shells  of  arbitrary 
shape.  In  the  middle  surface,  the  strains  are  given  by 

en  ■  v-[vMl + v2»2i  •  *! + M + i> 

*22  ’  cy^[aiU2,l2  +  C^,I1U11  '  Rj  +  2(4  +  1 

*12  ‘  2^V2,H  +  -  °iai,S2  -  »2a2.eiJ  +  2  *A 

where  the  rotations  ^  and  <J>2  are 

*  “  —  +  —  H 

"  Rj  +  Oj  ",s1 


and  the  rotation  about  the  norma.l  to  the  middle  surface  is 


\  2a1a2t(a2U2),i1  "  (cW»g2]  (3) 
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The  general  curvature-displacement  relations  used  by  Sanders  are  given 
by 

Kn  =  5^,0  +  Vi,$2J 

’*22-3^^,52“l  +  «i«!,i1]  w 

*12  *  +  '  *l“l.l2  -  ^,l1!  +  2  (r2  -  lj  % 

Equations  (1-4)  are  valid  under  the  following  assumptions: 

1.  The  Kirchoff-Love  hypothesis  holds;  that  is,  a  straight  line 
segment  that  is  perpendicular  to  the  middle  surface  of  the 
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deformed  shell  remains  perpendicular  to  the  deformed  middle 
surface  while  undergoing  negligible  strain  relative  to  its 
original  length 

2.  The  strains  and  rotations  of  line  elements  on  the  shell  mid¬ 
surface  induced  by  the  deformation  remain  small  in  comparison 
to  unity,  although  the  components  of  the  displacement  are  not 
necessarily  "small". 

Thus,  the  total  strains  in  terms  of  displacements  are 


C11  = 


11 


+  z  v. 


11 


-■*  sss 

22 
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22 


pk  = 
£12 


€12  +  22  *12 


By  choosing  the  coordinates  for  a  truncated  conical  shell  as  shown  in 
Figure  1,  the  general  parameters  appearing  in  equations  ( 1 ) - (4 )  may  be 


written 

as 

*1 

= 

X 

axial  coordinate  along 

a  generator 

^2 
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circumferential  coordinate 

U1 
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displacement  component 

in  x  direction 
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displacement  component 
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The  rotations  given  by  equations  (2)  and  (3)  consequently  reduce  to 


vcos  a  +  w. 


xsin  a 


«  -ifc+v 

vn  2 \x  ’x  xnln  a/ 


and  the  nonlinear  strain  and  curvature  relations  for  a  conical  shell 


become 
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(12) 


1  *  * 
w,  +  r  v,  co*  a  w,  +  v  cos  a 

m  xq>  2  *x _  _  tp  _ 

*Sqp  x  sin  a  2  .  _ 

“  x  sin  a 

A*  noted  in  the  introduction,  previous  authors  studying  conical  shell 
behavior,  notably  Seide^'11  and  Singer10,  have  employed  simplified  fornu 
lations  of  equations  (11)  and  (12)  by  neglecting  the  circumferential 
displacement  component  and  its  derivatives  in  the  curvature-displacement 
relations.  These  terms  were  also  neglected  in  the  nonlinear  strain- 
displacement  relations.* 


*  9 

Throughout  this  report,  the  Seide  formulation  will  be  referred  to  as 
"classical1'  conical  shell  theory,  and  the  terms  which  are  neglected  in 
said  formulation  will  be  traced  by  a  (*)  superscript,  as  in  equations 
(11)  and  (12). 
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TOTAL  POTENTIAL  FNERGY 


The  total  potential  energy  of  the  shell  in  the  presence  of  applied 

external  loads  may  be  expressed  as 


UT  = 


U  + 
m 


Ub  + 


V  + 
tn 


(13) 


where  U  is  the  strain  energy  due  to  middle -surface  stretching,  U,  is 
m  b 

the  strain  energy  due  to  bending,  is  the  potential  of  the  applied 

surface  loads,  V  is  the  potential  cf  the  applied  edge  loads,  and  V, 
m  b 

is  the  potential  of  the  applied  edge  moments  and  shears. 

An  analysis  of  the  nonlinear  behavior  of  a  conical  shell  in  the  presence 
of  external  loads  is  presented  in  Appendix  I;  in  Appendix  II,  the 
resulting  general  stability  equations  of  a  conical  shell  are  discussed. 

In  this  section,  however,  external  loads  are  not  considered;  instead,  the 
equations  of  motion  for  the  nonlinear  free  vibrations  of  a  conical  shell 
are  developed. 

In  the  absence  of  applied  external  loads,  the  total  potential  energy  of 
the  shell  is  given  by  the  strain  energy 

U  =  U  +  U,  (14) 

m  d 


The  strain  energy  due  to  middle-surface  stretching  may  be  expressed  as 


U 

m 


■  i/Ke* 


A 


+  N  e  +  N  7  )dA 

Cp  Cp  Xcp  Xcp 


(15) 


where 


N  ,  N 
x 


9 


and  N 


'  xcp 


are  the  middle- surface  forces  shown  in  the 


cp’ 


and  7  are  the  middle-surface  strains; 
xcp  * 


figure  on  page  11;  e^, 

and  A  is  the  surface  area. 

Likewise,  the  strain  energy  due  to  bending  may  be  expressed  as 


U 


b 


2M  h  )dA 
xcp  xcp 


(16) 


where  M  ,  M  ,  and  M  are  the  bending  and  twisting  moments  shown  in 
x’  cp'  Xcp 

the  figure  on  page  11,  and  x  ,  x  ,  and  x  3re  the  middle-surface  curvatures. 

x’  cp  xcp 
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KINETIC  ENERGY 


The  kinetic  energy  of  the  shell  is  the  sum  of  the  kinetic  energies 
associated  with  the  axial,  tangential  and  radial  velocities;  that  is, 

«PZ  *2 

T  *  j  ph  J  f  (u2  +  v2  -t  w2)x  sin  Ct  dx  dqj  (24) 

*1*1 

where  q^  and  cp2  taken  from  cpj  *  0  to  =  2  it  for  the  entire 

shell. 
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VARIATIONAL  PRINCIPLE 


The  application  of  Hamilton's  principle  requires  that  the  simultaneous 
first-order  change  in  the  Lagrangian  (£  =  T  -  U),  integrated  over  a 
specified  time  interval  with  respect  to  admissible  variations  in  the 
degrees  of  freedom  characterizing  the  state  of  strain  (namely,  u,  v,  and 
w),  must  vanish;  therefore. 


A- 


6u  |  (T  -  U)dt 
L1 
h 

6VJ  (T  -  U)dt 


/< 


6  /  (T  -  U)dt 

w 


=  o 


(25) 


or 


y 6ui  dt .  fiv  dt  -  o 


/ 6VT  it  -  o  > 


A1  dt  -  Au  dt  *  ° 


(26) 
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DIFFERENTIAL  EQUATIONS  OF  MOTION 


Appropriate  variation  of  the  kinetic  energy  [equation  (24)]  yields 


6  T  « 
u 

ph J/ufiu  x 

sin 

Ct  dxdcp  > 

q>  x 

6yT  " 

ph  J Jvfiv  x 

sin 

0!  dxdcp  > 

(p  X 

6  T  - 
w 

phjjwfiw  x 

sin 

Ci  dxd  cp  ^ 

cp  x 


whereas  integration  with  respect  to  tiae  gives 


fij dt 

t 

fc <** 

t 

At  dt 


Ph^iufiu]  2  x 

sin  a  dxd qj  - 

-///■ 

sin  Ct  U6u 

-\ 

dxdcpdt 

Cp  X  C1 

t  cp  x 

Ph  JJ  vfiv]  2  x 

sin  a  dxdcp  - 

sin  a  Vcpv 

dxdcpdt 

>(28) 

cp  x 

t  cp  X 

phj^fwfiw]  2  x 

sin  Ct  dxdcp  - 

Ph [ff* 

sin  a  ii6w 

dxdcpdt 

J 

cp  X  1 

t  cp  x 

The  variation  of  the  total  strain  energy  U  is 

6U  =  6U  +  (29) 

m  d 

where  6(  )  -  6  (  ).  Substitution  of  equations  (15)  and  (16)  into 

u,v,w 

equation  (29)  yields 


6U  =  ^  //[( N  6e  +  e  fiN  )+(N  fie  +  e  6N  )+(N  67  +  7  6«  )]dA 

2  JJ  xu  x  xu  x  cp  cp  cp  cp  xcp  *cp  xcp  xcp 


(continued) 
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-  b  /*/[(M  6h  +  H  6m  )+(rf  6k  +  k  6M  )  -  2(M  6h  ^  +  k  6M  > ]dA 
IJJ'-  xx  x  x  '  9  q>  9  9  x<p  x?  x9  xcp  J 

A 

(30) 

Now,  for  a  linearly  elastic  material  the  following  relations  hold  for  all 
force  and  moment  components. 

N6c  *  fc6N 

m6k  »  k6m  (31) 

Consequently,  equation  (30)  for  the  variation  of  the  strain  energy  may 
be  simplified  to 


6U  =// (N  6e  +  N Jk  +N  fir  -  M  6k  -  M  6k  +  2M  6k  )dA  (32) 
JJ  x  x  9  (p  xcp  xcp  xx  q>  9  xcp  xcp  '  ' 

A 

Substitution  of  equations  (11)  and  (12)  into  equation  (32),  with 
appropriate  integration  by  parts  and  rearrangement  of  terms,  yields 


6 

_ v 

sin 


U  fC\l  M  2M  \  /n 

t  a  JJ  toe  sin  a  x9,x  x  J  \sin  a 

to  x  *- 

*  * 

(vN  cos  a  +  w,  N  .  \ 

^"x'  sin'a  **  +  %w-,)cot  a  6v 


+  2N  +  xN 
xcp  X9 


J 


+  /[(xN  +  M  cot  a)6v]  2 

J L  xcp  X9 

9 


-  f\l  N  M  cot  a\ 

2  d  cp  +  /  M - 2- -  6v 

l  ^  J  \sxn  a  x  sin  a  J 

v 


9n 


«Pl 


dx 


(34) 
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5  L- 

W 

sin 


//lh- 


M 


cc.99 


x  sin  a 


2M  „  2M 

*».*9  .  h _ 5Sb$_ 

sin  a  cp,x  x  sin  Ct 


+  N  cot  .*  +  ~r —  [(w.  +  v  cos  a)N  +  x  sin  Ot  w,  N  1, 

V  sin  u  *■'  ’cp  xcp  ’x  xJ’x 


+  — 


x  sin  ^ 


* 


(w,  +  v  cos  Q)N  +  x  sin  Ot  «,  N  cp 

cp  cp  *x  XT 


/6w  dxdcp 


’9 


r\  2M  N 

f  !(xM  ),  -  M  -  — (v,  + 

J  |  j  x  x  cp  sin  a  sm  a  ’cp 


v  cos  a)  +  xw,  N  6w 
’x  x  | 


-  xM  6w, 

X  X 


dcp 


J  ||  sin  a  \x  si 


2M 

cp _ m  _ 

in  a  x 


2M  I 
xcp,x/ 


N 


£ 


— 1 , - (w,  +  v  cos  a) 

i„2  rv  9 


x  sin  a 


N  w,  I 
xcp  xl, 

sin  u  j 


+  -  -  V  --  6w, 


.  2 
x  sin  a 


9 


,CP’  92  x 

+  2[(M  )  Z]  Z 

cp.  1 


(35) 


EQUATIONS  OF  MOTION  IN  TERMS  OF  FORCES.  MOMENTS.  AMD  DISPLACEMENTS 

When  equations  (23) ,  (33),  (34),  and  (35)  are,  in  turn,  substituted  into  the 
variational  formulation  given  by  equations  (26),  three  nonlinear  partial 
differential  equations  of  motion  (Euler  equations)  are  obtained  for  the 
conical  shell  in  terms  of  forces,  moments,  and  displacements: 


N  -  N  N 

n  +  -£ - 2  +  -sua-  =  phU 

x,x  x  x  sin  a 
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(36) 


4 


I  W’X  1  I 

M  =  -D.v,  +  .  - +  — - - — (v,  +  v,  cos  a)  , 

/  xx  x  2  .2  ’cpcp  ’cp  | 

1  x  sm  ^  T  x  * 


!  1 


|y  ^  ^  | 

\  ■  -Drr +  r~r  (”-w  +  v>9  ~,s  “> +  vu'*x 

T  '  x  sin  a  ) 


1  *  * 
jw,  +  ^  V,  COS  Cl  W,  +  V  COS  a) 
1  ’xcp  2  ’x  _  ! 


M  =  D(l-v)  - ^r— 7- 

xcp  |  X  sin  U 


2  .  j 

x  sm  Ct  / 


Next,  the  three  equations  of  motion  given  by  equations  (36)-(38)  may  be 
written  in  terms  of  forces  and  displacements  by  virtue  of  equations  (40); 


that  is, 


N  -  N  N 

N  +  -* - 2  + -  phii  (41) 

x,x  x  x  sin  Cx 

*  *  * 

+  N  +  Sa  ,r .  ...P-g.ps-g  .  « 

x  sin  Ci  XV, *  x  X2  sln2  Q  [  '*x<p  x  ,2  si„2  0 


1  /  v,  \  /w,  +  V  cos  Cl\ 

if  V,  -  — -)  cos  a  -^-g  - I  N 

2  y  xx  f  x  \  x  sin  a  /  (p 


+  w,  N  =  phv 
’x  xcp  K 


„  „4  .  D  cos  Ci 

D  V  w  -r  — v,  x 

c  2.2. 
x  sm  u 


ic  Jc 

(V2vJ,  -4  (-^)  ] - t— rjN 

c  J’cp  \  x  /,  I  x  sin  a  j  cp 


cos  OL 


-r(w,  +  v  cos  Ct)N  +  x  sin  OL  w,  N  ], 

xcp  ’x  xJ’x 


W,  .  +  V  COS  Cl  \  I 

— ^ — — — -  ]  N  +  w,  N  I  =  -phw 

x  sm  Ci  /  cp  x  xcpl,  ) 


where  the  harmonic  ooeratcr  V  is  defined  as 

c 


V  }  (  );xx  +  x(  }’x  +  2  .2  (  }  * cpcp 

x  fi\n  CL  ~~ 
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EQUATIONS  OF  MOTION  IN  TERMS  OF  DISPLACEMENTS 

Equations  (41;-(43)  may  further  be  written  in  terms  of  displacements  only 
by  utilizing  equations  (39).  The  displacement  equations  of  motion  thus 


become 


JlJ|  xi  .  ,  %  ~  “  cos  g ,  i  2  : 

,  2i|u*x  2  W*x  V  x  x  sin  Ci  -2.2  ^w*cp  v 

i-v  (<  2x  sm  a  Y 


* 

2  ~  2 
ccs  a 


*  II  i_v  {  1  o  v,  -  w  cos  a 

+  2vw,  cos  OL  j  H - j  u,  +  —  w, - - £P — - — — - 

’cp  J),x  x  (  ’x  2  ’x  x  x  sin  a 


*  -k  \ 

1  2  2  ~  2  / 
— r - 2 - (w,  +  v  cos  a  +  2vw,  cos  a)) 

2x/  sin  a  Y  ^  ) 


2x  sin  a  |  v’x  +  x  sin  a  x  T  x  sin  alw,cp 


v  .  ’x  , 

-  7  +  +  v  cos  a/ 


-  pii  (46) 


_  -  [v,  -  W  COS  a  , 

-JL) _ 1 _ Is _ +  E+ _ L_ 

,  2  x  sin  Ct  x  sin  a  x  „  2 

1-v  (  2x  si 


2  2  ~  2 

- = - (w,  +  v  cos  a 

.  2  _  cp 
sin  O. 


*  i  2  i_i  u» 

+  2vw,  cos  a)  +  v(u,  +  -z  w,  )  4-  v,  4-  — — -  -  — 

’cp  ’x  2  ’x  J  ,(p  2  ’x  x  sin  a  x 


w,  *  .  u,  w, 

4 - 7 — “(w,  4  v  cos  Ct)  4-  — -  v,  4 - -  -  —  4 - ■ 

x  sin  CT  ’cp  x  ’x  x  sin  a  x  x  sin  a 


h  cos  a 


(%+VCOS  a)|+“r“2“W’xxcp+  X 


12  sin  a 


*  *  * 
w,  w,  4-  v.  cos  a 

'  vm  '  ecrntci  •  rrvr\ 


2  2 
x  sin  a 


1-v.  v’x. 

4-  — r~(v, - )cos 

2  xx  x 


J  .  ct  a|(“-cp+  v  c°s  g\  [»  VV 
x  |  \  x  sin  a  /  v  x  si 


-  wcosCt 


x  sin  a 
(continued) 
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+  — 2~~~2~(w4  +  v2c°s2a  +  2vw’c  cos  a  cos  a)+  v^u»x  +  2W*x) 

2x  sin  u  J 

w,  (l-i  )  u,  w,  *  ) 

x  cp  v  .  x  ,  1 

+  - z -  v,  + - .  —  -  —  H - : - r(w,  +  v  cos  ci)  i  =  pv 

2  xx  sLn  u  x  x  sin  u  ’cp  I 


_4  D  cos  u  ,„2  ^ 

7  w  +  -r - - —  (V  v) , 

c  2.2  c  ’cp 

x  sin  u  T 


..  ■  >fi 


x  ’x*  (l-v2)x  sin  a 


cos  G£; 


V,  -  W  cos  a  2  2*2  * 

+  — * — : —  - +  — r - z — (w,  +  v  cos  a  +  2vw,  cos  G) 

x  sin  Q  _  2  .  2  '  cp  *tp 

9v  Cin  O  “  “ 


2x  sin  Ci 


12)1  *  1-v  u,cd 

+  v(u,  +  —  w,  )  +  (v.  +  v  cos  a)  — t—  v,  +  — -r-— - 

x  2  ’x  |  i  ’cp  2  ’x  x  sin  a  x 

W>x  *  I  1  2 

+  - :  - . : (w,  +  v  cos  a)  +  x  sin  G  w,  ;  u,  +  —  w, 

x  sm  a'  ’cp  *x  ’x  2  ’x 


v .  -  w  cos  a 

cp  u 

+  v  — * — : — - +  -  + 


:<  sin  u 


*  J 

+  2 vw,  cos  Ci) 

’cp 


u  .  1  ,  2  .  2  ~  2  ,, 

x+  ,"r  .  2  ;(”■<, +  v  cos  a 

2x  sin  a  T 


|  +  l/v’cp+  v  c°s  a\[v»cp  ~  w 

}  ’x  | \  x  sin  a  /  x  sin 


cos  a 


u  1  2  2  ~  2  * 

+  —  +  — z - - — Cw,  +  v  CDS  a  +  2vw,  cos  Q) 

x  ,,  2  .  2  „  ’cp  ’cp 

2x  sm  Ci  T  Y 


w,  (1-v) 


+  V(~U>y.  +  2  W»X)  + - 2" 


u,  w, 

V  +  - JBL-  _  X  +  - !*_(w 

’x  x  sin  Cl  x  x  sin  Cl''  ’cp 


+  v  cos  Cc)  /  =  -  phw  (48) 
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BOUNDARY  CONDITIONS 

The  associated  boundary  conditions  which  result  from  application  of  the 
variational  principle  [equations  (25)]  are  given  in  terms  of  forces, 
moments,  and  displacements  as  follows: 

(a)  along  x  =  x^  and  x  =  x^ 


or  6w  =  0  (51) 


These  boundary  conditions  may  be  written  also  in  terms  of  forces  and 
displacements  by  using  equations  (40);  thus. 
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;a)  along  x  =  x^  and  x  = 


N  =  C 
x 


or  6u  *  0 


.  D(  1- .  )  cos  ~i,  *  .  1 

■i  +  — ~ — —z - (w,  +  -z  v,  cos  Q  - 

xc  2.2.  xcp  2  ’x 

x  sm  G  T 


*  * 
w,  +  v  cos  G 

-* - )  =  o 

x 

or  6v  »  0 


i  j  w’x  1  *  111 

Y sln  •'*  i  "•*x + '  ~ +  7Y77("’w +  v>9 cos  a>  i 


+D  sin  -~t -  +  — - z - (w,  +  vs  cos  G)  +  v  w,  ] 

'x  2  .  2  „  *cpcp  °©  ’xxJ 

x  sin  Lt  ^  T 


20(1-; )  ^  1  .  "■»  T  CM  a 
x  sin  G  w’xqp  ^  V’x  C0S  x 


+  x  sin  Ct  w  N  +  N  (w,  +  v  cos  a)  *  0  or  6w  *  0 

’x  x  xqr  *q>  — 


*1 

w,  +  v\ — -  +  — - = — (•»,  _  +  v,  cos  a)]  =  o 

xX  ‘  x  X2  Sin2  a  &  * 


or  Gw,x  *  0 


(b)  along  cp  *  ^  and  cp  =  cp2 


N  =0 
xcp 


or  6u  =  0 


,  D  cot  u  ’x  .  1  *  * 

N  +  -  - +  —Z - z — IW,  +  v,  cos  Q) 

w  x  x  2,2.  ’cpcp  cp 

x  sin  a  ^  T 


vw,  =  0 
’xx 


or  6v  =  0 


w 

•“,*  +  %  cos  a>  +  ‘'“■xx  , 

x  sm  u  T  ’cp 


-  2D( l-v )  w,  +  £  v,  cos  a 
xcp  2  x 


W,  +  V  COS  a 


*x  (continued) 


(53) 


REDUCTION  TO  CLASSICAL  EQUATIONS 


A  simplified  set  of  equations  of  motion  may  be  obtained  if  simplified 

strain  and  curvature  relations  similar  to  the  well-known  von  Karman- 

Donnell  relations  for  circular  cylindrical  shells  are  employed  instead  of 

the  more  accurate  relations  deduced  in  this  study  [equations  (li)  and 

9 

(12)  ].  These  simplified  relations  are  the  same  as  those  used  oy  Seide 

and  Singer*^  for  the  stability  analysis  of  conical  shells;  they  may  be 

* 

obtained  by  neglecting  all  terms  marked  by  a  (  )  superscript  in  equations 
(11)  and  (12).  Thus, 


1  2 

c  =  U,  +  —  W, 

x  x  2  x 


w  cos  O. 


€  =  -  + 
cp  X 


x  sin  Q 


2  2 
2x  sin  a 


y  =V  +  - JB- 

xcp  *x  x  sin  a  x  x  sin  a 


k  =  w, 

X  XX 


Kcp  2  2  x 

Y  x  sin  O. 


xcp  x  sin  0:  2  . 

T  x  sin  0£  J 

The  force-displacement  nonlinear  equations  of  motion  resulting  from  sub¬ 
stitution  of  equations  (55)  and  (56)  into  the  Lagrangian  and  carrying  out 
the  variations  in  accordance  with  Hamilton's  principle  may  be  obtained 

•k 

here  simply  by  neglecting  the  terms  marked  with  a  (  )  superscript  in 


equations  ( 36 ) - ( 38 ) ;  that  is, 


N  -  N  N 

N  +  -* - 2  +  -  phU 

x,x  x  x  sin  Ct 


N  2N 

+  n  +  -SB  =  Phv 
x  sin  Q  xcp,x  x 


(continued) 
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(57) 


Equations  (58)  agree  vith  those  presented  in  Reference  25  for  a  nonlinear 
postbuckling  analysis  of  a  conical  shell. 

The  three  equations  of  motion  given  by  equations  (58)  may  be  reduced  in 
number  by  one  through  introduction  of  a  stress  function  F  defined  by 
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N  F, 

x  x 


x  2  2 

x  sin  Q 


N 

-se  -  F 

u  1  ». 


N  ,  F, 

-22  =  - I - (F  .  -12) 

h  x  sin  Q  *xcp  x 

such  that  the  first  two  equations  of  (58)  are  satisfied  identically. 

To  establish  a  relationship  between  F  and  w  independent  of  u  and  v, 
the  strain  compatibility  equation  for  a  conical  shell,  as  employed,  for 

q 

example,  by  Seide  , 

r”P.f e  .  ifs6*  +  ix£  .  r<w)  (60) 

x2  sin  Q  *  sln  °  X2  sin2  Q  ’’•xx  * 

is  used  (where  T  is  a  function  of  w).  Substitution  of  equations  (55)  into 
equation  (69)  gives 

2 

w,  cos  a  /w>  \  9 

F(w)  ~  — — — —  -  — ~ - = —  w,  -  w,  w  (  — x  sin  a  w,  w,  , 

x  si.i  Q  x2  a  *xcp  ’xx  cpcp\  x  /  *x  ’x^ 

(60a) 

Substitution  of  the  constitutive  law  [equation  (17)]  and  equation  (59) 
into  equation  (60),  in  conjunction  with  equation  (60a),  then  yields  the 
desired  compatibility  equation  in  terms  of  the  stress  function  and  radial 
displacement  only  as 


E 

V  F  =  —z - z — 

c  2.2. 
x  sin  Cc 


2  rSg\  .  2 

w,  -  w,  w,  -l — xl  -  x  sin  Q  w,  w, 

xcp  XX  CfCp  \  X  /,  X  XX 


Ev,  cos  C* 

XX _ 

x  sin  a 


The  lateral  equilibrium  equation  may  be  written  in  terms  of  the  stress 
function  F  and  the  radial  displacement  w  by  substitution  of  equations 
(59)  into  the  third  of  equations  (58);  thus. 
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Equations  (61)  and  (62)  together  represent  the  s imp lifted  nonlinear 
vibrational  behavior  of  a  conical  shell.  They  are  the  conical  shell 
counterparts  of  the  von  Karman -Donne 11  equations  for  circular  cylindrical 
shells.  When  equations  (61)  and  (62)  are  linearized,  the  resulting 
equations  are  those  describing  the  classical  behavior  of  a  conical  shell 
in  free  vibration;  thus. 


D74w  -  — -?t . °  F,  +  phw 

C  X  X* 


and 


Ew,  cot  a 
xx _ 

X 


(64) 


These  linear  equations  were  solved  by  Seide  in  Reference  11. 
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SUGGESTED  METHODS  OF  SOLUTION 


The  differential  equations  (46) ,  (47) ,  and  (48)  reflect  the  use  of  the 
complete  strain  and  curvature  relations  given  by  equations  (11)  and  (1?.); 
they  are  three  highly  nonlinear  coupled  equations  in  the  three  dependent 
variables  u,  v,  and  w.  On  the  other  hand,  the  simplified  "classical" 
strain  and  curvature  relations  given  by  equations  (55) ,  with  the  midsur¬ 
face  inertia  terms  neglected,  admit  reduction  to  two  simultaneous 
differential  equations  involving  only  the  stress  function  F  and  the 
radial  displacement  w.  This  reduced  set  of  two  simultaneous  equations 
[equations  (61)and(62)]  is  still  highly  nonlinear  and  coupled  but  probably 
more  manageable  in  the  classical  sense  than  the  initial  set  of  three 
given  by  equations  (41)-(43).  However,  the  classical  equations  are  valid 
only  for  motion  of  the  shell  in  which  a  large  number  of  circumferential 
waves  is  present. 

As  a  first  step  toward  solving  the  conical  shell  equations  described 

above,  it  would  be  feasible  to  employ  properly  modified  techniques  used 

previously  for  the  limiting  case  of  a  circular  cylindrical  shell  (as 

18 

discussed  by  Mayers  and  Wrenn  ) . 

All  of  the  existing  analytical  studies  of  the  nonlinear  vibrations  of 
thin  circular  cylindrical  shells,  except  for  those  reported  in  References 
18,  21,  and  22,  have  the  inherent  assumptions  (1)  that  the  radial  displace¬ 
ment  is  representable  by  a  chessboard  pattern,  (2)  that  the  shell  possesses 
many  circumferential  waves  (i.e.,  shallow-shell  behavior),  and  (3)  that 
the  midsurface  inertia  terms  can  be  neglected.  In  accordance  with 
assumptions  (2)  and  (3) ,  all  studies  have  involved  the  solution  of  the 
classical  equations  [equations  (103)  of  Appendix  III]  based  upon  von 
Karman-Donnell  theory*.  Assumption  (1),  however,  seems  to  be  intuitively 
incorrect.  As  noted  in  References  18  and  21,  use  of  the  chessboard 
pattern  for  the  radial  displacement  leads  to  a  violation  of  the  periodicity 


-'Mayers  and  Wrennl8  also  carried  out  a  first-order  approximation  solution 
of  the  equation  set  [equation  (98)]  based  upon  Sander's  strain  terms 
with  midsurface  inertia  terms  retained. 
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condition  for  the  circumferential  displacement  v  given  by 


or. 


v(x,cp,t)  =  v(x,q),2«  ,t) 


(65) 


2n 

j'v,cpdcp=0  (66) 

0 

when  finite  displacements  occur. 

The  fixed -parameter  solution  for  a  circular  cylindrical  shell,  first 

22 

presented  by  Evensen  on  the  basis  of  the  Galerkin  procedure  and  again 
by  Mayers  and  Wrenn^  as  a  special  case  using  the  Rayleigh-Ritz  approach, 
may  be  extended  to  the  conical  shell  case  by  assuming  a  radial  displace¬ 
ment  which  satisfies  both  the  geometric  boundary  conditions  and  the 
periodicity  condition. 

If  the  geometric  boundary  conditions  are  taken  to  be  those  for  classical 
simple  supports,  namely 


w(x,cp,t)  ■  w(x2 ,  cp,  t )  =  0 


v(x,cp,t)  =  v(x2,cp,t)  =  0, 


(67) 


then  a  corresponding  second-order  approximation fbr  the  radial  displace¬ 
ment  w  which  satisfies  the  boundary  and  periodicity  conditions  may  be 
taken  as 


.  .  .  itmCx-x.) 

w  =  A^(t)  sin  _ 1_  cos  n(p  + 


x„-x. 

i.  1 


2.2.  . 
n  AL(t) 

4xsina  cosCt 


.  2  mrt(x-x,) 
sin  1 


VX1 


(68) 


where  A^(t)  is  an  arbitrary  function  of  time.  However,  by  restricting 
the  time-dependent  coefficients  to  terms  proportional  to  A^(t)  and 
2 

A^(t),  the  resultant  motion  is  forced  to  be  periodic,  although  not 
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necessarily  harmonic. 

A  Ray le igh-Ritz  approach  (as  described  in  Reference  18  for  a  circular 
cylir.  liical  shell)  may  then  be  used  by  substitution  of  the  assumed 
deflected  shape  into  the  compatibility  equation  [equation  (61)]  and 
obtaining  a  particular  solution  for  the  stress  function  F.  The  stress 
function  thus  obtained  and  the  assumed  deflection  function  w  should  lead 
to  satisfaction  of  the  requirement  that  the  circumferential  displacement 
v  be  periodic  in  2it.  F  and  w  are  then  substituted  into  the  Lagrangian 
written  in  terms  of  stresses  (and  therefore  in  terms  of  the  stress 
function)  for  midsurface  behavior  and  in  terms  of  w  for  bending  behavior. 
The  application  of  the  variational  principle  relative  to  Aj(t)  would 
thus  yield  a  second-order  nonlinear  differential  equation,  which  through 
a  single  change  of  variable  may  be  reduced  to  two  coupled  first-order 
nonlinear  differential  equations.  These  may  then  be  solved  by  a  standard 
Runga-Kut.ra  numerical  technique. 

The  solution  obtained  in  the  preceding  manner  will,  as  previously  stated, 

18 

be  periodic.  Mayers  and  Wrenn  waived  this  restriction  for  a  circular 

cylindrical  shell  by  allowing  for  a  free-parameter  solution.  The 

22 

boundary  conditions  were  the  same  as  those  stipulated  by  Evensen 
(free  support).  The  boundary  conditions  will  not  significantly  influence 
the  frequencies  of  vibration  when  a  large  number  of  waves  are  present  in 
the  axial  direction;  thus,  for  long  shells,  on  a  minimum  energy  basis, 
the  nonperiodic  behavior  governs. 

Extending  the  approach  to  a  conical  shell,  an  assumed  radial  deflection 
shape  may  be  taken  as 

nm(x-x^)  2rot(x-Xj) 

w  =  A  (t)  sin  -  cos  mp  +  A_(t)  cos  -  +  A  (t)  (69) 

1  x2~*l  ^  X2*x^  j 

where  A^(t),  A^(t)  and  A^(t)  are  arbitrary  functions  of  time. 

A  displacement  function  of  this  type  was  also  used  in  a  study  of  the 
postbuckling  behavior  of  conical  shells  presented  in  Reference  25. 

Next,  using  the  assumed  displacement  function  w,  the  solution  for  the 
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stress  function  F  is  obtained  from  the  compatibility  equation 
[equation  (61)],  and  the  parameter  A^(t)  is  then  determined  by  impos¬ 
ing  the  continuity  condition  on  the  circumferential  displacement  v.  The 
v  and  F  functions  thus  obtained  are  finally  substituted  into  the 
Lagrangian  (in  a  manner  similar  to  that  described  previously  for  the 
fixed-parameter  case)  and  variations  are  taken  with  respect  to  A^(t) 
and  A^u).  The  resulting  nonlinear  differential  equations  may  again 
be  solved  using  a  Runga-Kutta  numerical  technique. 

Application  of  this  method  to  the  circular  cylindrical  shell  resulted 
in  a  solution  which  displayed  nonperiodic  motion;  the  same  type  of 
behavior  can  also  be  expected  in  the  case  of  a  long  conical  shell. 

The  two  solution  approaches  described  above  involve  much  greater  diffi¬ 
culty  vhen  applied  to  a  conical  shell  than  when  applied  to  a  circular 
cylinder,  since (1)  the  classical  equations  of  equilibrium  [equations  (61) 
and  (62) ]  for  a  conical  shell  are  more  complicated  than  those  for  a 
circular  cylindrical  shell  [equations  (106)  and  (107)  of  Appendix  III], 
and  (2)  the  stress-compatibility  equation  for  the  conical  shell  is  a 
bihansonic  equation  with  variable  coefficients.  Nonetheless,  the  methods 
of  solution  may  still  be  used  and  they  are  presently  undergoing  further 
study  in  connection  with  the  overall  research  effort  related  to  conical 
shells. 

When  the  number  of  circumferential  waves  is  not  large,  the  classical 
equations  are  invalid  and  should  not  be  used;  rather,  the  equations  of 
motion  [equations  (4l)-(43)]  based  on  the  general  strains  and  curvatures 
[equations  (11)  and  (12)]  with  retained  midsurface  inertias  should  be 
solved.  Alternatively,  a  direct  approach  may  be  used  by  substituting 
assumed  time-dependent  displacement  functions  directly  into  the  total 
potential  energy  as  given  by  equation  (23)  and  applying  the  variational 
principle  relative  to  the  assumed  time-dependent  parameters. 

18 

Mayers  and  Wrenn  used  the  latter  procedure  for  a  circular  cylindrical 
shell;  however,  due  to  the  obvious  algebraic  difficulties,  only  one 
term  in  the  radial  displacement  function  was  retained,  thus  eliminating 
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the  possibility  of  the  shell  deforming  in  its  typically  nonlinear  fashion. 
As  a  result  of  this  forced  chessboard  pattern,  the  resulting  motion  was 
periodic;  however,  for  very  small  amplitudes  of  vib1  •'.tion,  the  periodic 
solution  is  an  excellent  approximation. 


When  extending  this  approach  to  the  case  of  a  conical  shell,  displacement 
functions  for  u,  v,  and  w  which  satisfy  the  geometric  boundary  conditions 
for  w  and  the  periodicity  condition  may  be  taken  as 


mrt  (x-x^) 

w  -  A, (t)  sin  - 

1  x2-xx 


cos  nt: 


A 


mrt(x-x^)  2nm(x-x^) 

u  =  A  (t>  cos  -  cos  nep  +  A  (t)  cos - cos  2ncp 

X2*x^  j  x^  ™x^ 

2nnt  (x-x  ) 

+  A. (t)  cos  - 

4'  x„-x. 


mu (x-x  )  2nnt(x-x^) 

v  =  A  (t)  sin  -  sin  nCp  +  A,(t)  sin  -  sin  2ncp 

D  X^-X,  D  X„ -X, 


1 

y  v70) 


2  *1 


2  “1 
+  A^(t)  sin  2ncp 


Substitution  of  equations  (70)  into  equations  (23)and(24)  and  the  enforce¬ 
ment  of  the  simultaneous  vanishing  of  the  first  variation  of  the 
Lagrangian  with  respect  to  the  various  A^(t)  functions  would  lead  to  a 
set  of  seven  second-order  coupled  nonlinear  differential  equations.  After 
a  change  of  the  independent  variable  and  subsequent  reduction  of  the 
seven  equations  to  fourteen  first-order  nonlinear  differential  equations, 
the  solution  is  then  obtained  using  a  standard  Kutta-Merson  numerical 
technique . 


The  solution  techniques  described  above  all  suffer  from  a  common 
disadvantage;  namely,  the  number  of  terms  which  may  be  retained  in 
the  assumed  displacement  functions  are  limited  as  a  result  of  algebraic 
complications.  However,  a  numerical  approach  is  presently  being  studied 
in  which  a  general  series  solution  for  the  various  displacements  may  be 
assumed  and  a  large  number  of  terms  retained.  The  method  may  best  be 
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understood  by  demonstrating  it  for  the  problem  of  the  nonlinear  vibra¬ 
tion  of  a  beam.  This  is  accomplished  in  Appendix  V. 

The  proposed  method  basically  involves  the  selection  of  assumed  displace 
ment  functions  in  the  general  series  form 


j l 

W  =  53  $ •  (x.O  A.  (t) 
i=l 

m 

V  =£  ft(x.cp)  A  (t)  (71) 

i=i+l 

n 

U  =  53  ^(x.cp)  ^(t) 
i-nrt-1 

where  A.(t),  i=l,2,...n  are  arbitrary  functions  of  time. 

Substitution  of  equations  (71)  into  equations  (23)  and  (24)  and  the 
enforcement  of  the  simultaneous  vanishing  of  the  first  variation  of  the 
Lagrangian  with  respect  to  the  various  A^(t)  (without  multiplying  out 
the  series  or  carrying  out  the  integrations)  will  lead  to  a  set  of  n 
second-order  coupled  nonlinear  differential  equations.  After  a  change 
of  the  independent  variable  and  subsequent  reduction  of  the  n  second- 
order  differential  equations,  the  solutions  mav  be  obtained  using  the 
numerical  procedure  demonstrated  in  Appendix  V. 

The  most  general  set  of  appropriate  displacements  for  a  conical  shell 
which  satisfies  the  simply  supported  conditions  and  periodicity  of  the 
circumferential  displacement  would  be  equations  (71)  with 

V 

<5.  (x,cp)  =  sin  ip(x-x  )  cos  jcp  ,  l£i £,& 

1  j=l 

m' 

<£.  (x,cp)  =  cos(i-£)(3  (x-x  )  cos  jcp  ,  i+l^rSm  (72) 

1  1  j-A'+l 

n ' 

4>.  (x,cp)  =  sin(i-m)P (x-x. )  sin  jcp  ,  mfl^i^n 

1  1  j=m'+l 


35 


vlu-  re 


P  - 


(73) 


la  osifa',  may  be  taken  as  any  type  of  function;  the  numerical 

technique  Joes  not  differentiate  between  types  of  functions  and  is  there¬ 
fore  not  limited  in  this  respect. 

The  main  merit  of  this  proposed  numerical  technique  is  that  there  is 
no  necessity  of  expanding,  multiplying,  and  integrating  the  assumed  dis¬ 
placement  functions  analytically  at  any  time  during  the  solution.  Further 
merits  and  some  drawbacks  of  the  method  are  discussed  in  Appendix  V. 

Finally,  another  method  of  solution,  which  involves  the  use  of  a 
mod  if ied-Reissner  variational  principle,  should  be  investigated. 

»his  method  was  introduced  in  Reference  28  and  again  used  successfully 
in  Reference  29.  It  involves  writing  the  strain  energy  of  the  system  in 
such  a  manner,  that  Hooke's  law  [equations  (17)]  is  not  necessarily 
enforced  prior  to  variation  for  middle-surface  stresses  and  strains 
though  it  is  forced  for  the  bending  terms.  This  approach  is  shown  to  be 
highly  convergent  and  accurate  for  both  thin  plate  and  shell  problems 
undergoing  finite  deflections  even  in  the  presence  of  inelastic  deforma- 
t ions . 


This  method  of  solution  and  all  the  others  previously  described  are 
presently  being  investigated.  The  results  and  conclusions  are  expected 
to  provide  the  same  level  of  knowledge  for  the  nonlinear  vibrations 
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CONCLL'S  ION'S 


The  governing  differential  equations  and  associated  boundary  conditions 
for  the  nonlinear  vibrations  of  a  thin,  truncated  conical  shell  have  been 
developed  by  means  of  Hamilton’s  principle.  The  equations  are  based  upon 
the  nonlinear  strain-displacement  relations  deduced  from  those  developed 
by  Sanders  for  thin  shells  of  arbitrary  shape. 

Governing  differential  equations  and  boundary  conditions  have  also  been 
developed  for  the  nonlinear  behavior  of  a  conical  shell  in  the  presence 
of  static  and  dynamic  loads.  These  equations  apply  equally  well  to 
forced  nonlinear  vibration  problems,  nonlinear  vibrations  of  a  prestressed 
shell  and  nonlinear  prebuckling  and  postbuckling  behavior  of  an  arbitrarily 
loaded  shell.  Linearized  stability  equations  and  boundary  conditions 
using  the  nonlinear  terms  of  this  study  have  also  been  developed. 

All  equations  derived  in  this  report  have  been  shown  to  reduce  (by  a  set 
of  simple  substitutions)  to  corresponding  equations  for  circular  cylin¬ 
drical  shells  and  circular  flat  plates.  The  equations  have  also  been 
shown  to  reduce  to  the  so-called  "classical"  conical  shell  equations 
when  the  circumferential  displacement  terms  in  the  curvature  and  nonlinear 
strain  relations  are  neglected  and  the  midsuiface  inertias  are  omitted. 

Various  methods  of  solution  extending  from  those  applied  in  similar 
investigations  of  circular  cylindrical  shells  have  been  presented,  and 
a  raodified-Reissner  variational  principle  me thod  of  solution  has  been 
suggested.  Finally,  a  new  direct  numerical  solution  scheme  has  been 
proposed  and  demonstrated  analytically. 

Future  efforcs  in  this  area  should  include  implementation  of  the  proposed 
methods  of  solution  for  the  nonlinear  vibrations  of  a  conical  shell  and 
a  comparison  of  the  various  methods.  The  modif ied-Reissner  variational 
principle  approach  and  the  new  numerical  scheme  should  at  first  be 
employed  to  obtain  solutions  for  the  nonlinear  vibrations  of  a  circular 
cylindrical  shell.  Comparison  with  the  previous  solutions  should 
determine  the  merits  (if  any)  of  these  methods,  and  it  may  then  be 
decided  whether  they  should  be  applied  to  the  conical  shell  problem. 
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Finally,  these  same  methods  of  solution  should  be  applied  to  nonlinear 
forced  vibration  problems  and  to  the  important  design  problems  of  the 
nonlinear  post buck ling  behavior  and  maximum  strength  of  initially 
Imperfect  conical  shells. 
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APPENDIX  I 


EQUATIONS  FOR  THE  NONLINEAR  BEHAVIOR  OF  CONICAL  SHELLS 
UNDER  ARBITRARY  LOADS 

Equations  for  the  nonlinear  behavior  of  conical  shells  acted  upon  by 
externally  applied  arbitrary  loads  may  be  derived  without  significant 
further  difficulty  by  adding  to  the  strain  energy  U  [equation  (14)] 
the  potential  of  the  applied  loads;  that  is, 

T  L  m  b  (74) 

where  U  is  the  total  potential  energy  of  the  shell  in  the  presence  of 
applied  loads.  V  is  the  potential  of  the  applied  surface  loads 

ij 

(qx,qr.  and  q^) ,  is  the  potential  of  the  applied  edge  loads  (N  ,  N^, 


m 


and  ^x(p)»  and  is  the  potential  of  the  applied  edge  moments  and 

shear  forces  (M  ,  M  ,  M  ,  Q  ,  and  Q  ) . 

x  <p  x<p  x  cp 

The  applied-load  potentials  may  be  represented  as 


VL  *  -JJ (q^u  +  q^v  +  qzw)  x  sinQ 


dxdcp 


cp  x 


(75) 


m 


x2  r  _  92 

=  -J  [(Nxu  +  Nx(pv)  x  sinO  ]  d9  -J  |y  +  Nycpu 

9  xj  x  9i 


dx 


(76) 


M  w, 


V  =-f[x  Sina  (-M  w,  +  X(P  +  Q  w) 

b  J  l  x  ’x  x  sin  a  x 


9 


-f  +m  w,  +  Q  w 

J  [  x  sm  Q  X9  x  9 


9, 


dx 


9j 


d9 


(77) 


The  required  variations  are 
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6  v  - 

u  m 


6  V  =  -  f  f  q  x  sinu  6u  dxdcp 
u  L  J  J  x 

V  x 

=  -  £  J*  q  x  sinu  6v  dxd9  ) 

c  X 

Vl  -  -//  q_^x  sinu  6w  dxd9  J 

9  x 

/■  _  xo  r-  % 

[xN  6uf  dcp+  ![N  fiul  dx 

v  J  '  rn 


1  X 


6  V  =  -  sin 
v  m 


6  V  =  0 

W  m 


/—  X2  /*  —  c? 2 

[xN  .  6v]  dtp  +  /  [N  ^  6v]  dx 

x.  J  ^  <•« 


1  x 


6  V  =  0 
u  b 


6  v,  =  o 

v  b 


6  V  *  sinu 
w  b 


X2 

/|xfer'5‘)S“|  d¥+/'(V*‘V 


i 

sinu  /  |xM.  6w,  dc  +  f |  M;p  6w, 

y  I  /  irfe  * 


-  .  %  *1 
dx  -  2[[M  6w]2  ]2 

CPi  Xj 


X1  *P1  T1 

Equations  (78)  through  (80)  are  added  into  the  variational  equations 
equations  (26)  ]  such  that 

6  /  (T  -  U  )  dt  =  6  f  (T  -  U  -  V  -  V  -  V  )  lit  =  u 

u.v.w7  r  u.v.w J  K  u  VL  m  U 


where  the  expressions  for  6  T  and  6  U  are  given  by  equations 

U,v,w  u ,  v ,  w  ^ 

(27.  and  (33) - (35 ) ,  respectively. 


The  variation  results  in  three  Euler  equations,  which  are  the  same  as 
those  given  by  equations  (36)-(38)  with  the  addition  of  the  surface  loads 
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q^,  q  and  q^,  respectively,  to  the  left-hand  side  of  each  equation. 

The  equations  of  motion  in  terms  of  stresses  and  displacements  are  the 
same  as  those  given  by  equations  (41)- (43)  with  the  respective  surface 
loads  added  to  the  left-hand  sides  of  the  first  two  equations  and  to  the 
right-hand  side  of  the  third;  the  same  is  true  for  the  displacement 
equations  of  motion  given  by  equations  (46)-(48). 

The  resulting  boundary  conditions  are 


(a)  along  x  =  Xj^  and  x  *  x^ 


N  =  N 
x  x 


M  cos  Ci 

XCp 


xcp  x  sin  Ci 


(xM  )  ,  -  M  - 

x  ’x  9 


sin  O: 


=  xQ 

x  sin  a 


M  =  M 
x  x 


or  6u  =  0 


or  6v  =  0 


+  — 7— —  (w,  +  v  cos  a)  +  xw,  N 

sin  a  *cp  xx 


or  6w  =  0 


or  6w,  =  0 

—  x 


(b)  along  9=cp^  and  9  =  9^ 


N  =  N 
X9  X9 


or  6u  =  0 


M  cos  a 


N - ; -  =  N 

9  x  sin  a  9 


or  6v  =  0 


2 (xM  ) , 
X9  x 


x  sin  a 


+  — — 7  (w,  +  v  cos  a) 

x  sin  u  9 


+  N  w.  =0  -  M 


or  =  n 


ffl 


<  c  * 


ac  the  corners  of  the  segment 


(84) 


M  =  M  or  fiw  1  0 

xr  xc  — 

All  equations  and  boundary  conditions  presented  in  this  appendix  reduce 

o  25 

to  those  of  classical  conical  shell  theory  '  by  linearization  of  the 
equations  and  neglect  of  those  terms  marked  by  a  (*)  superscript. 
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APPENDIX  II 


GOVERNING  EQUATIONS  AND  BOUNDARY  CONDITIONS  FOR  THE 
STABILITY  OF  CONICAL  SHELLS 

Additional  definition  of  symbols: 


M  , 

X 

M  , 
9 

** 

-c 

additional  moments  causeo  by  buckling 

M 

,  M 

,  M 

moments  prior  to  buckling 

X 

o 

% 

N  , 

X 

N  , 

9 

N 

xcp 

additional  f-  oes  caused  by  buckling 

N 

,  N 

•  NV 

forces  prior  to  buckling 

X 

o 

,  x9o 

e  , 

e  , 

y 

linear  strain  terms 

X 

9 

xcp 

e' , 

e ' , 

y' 

nonlinear  strain  terms 

X 

9 

xcp 

Linearized  stability  equations  and  boundary  conditions  for  a  conical 
shell  may  be  readily  obtained  from  the  general  equations  of  nonlinear 
behavior  of  conical  shells  derived  in  Appendix  I. 


The  equilibrium  at  the  time  of  buckling  is  obtained  by  consideration  of 
the  additional  work  done  during  buckling.  The  displacements  u,  v,  and 
w  are  now  the  additional  dis  Pi  acements  caused  by  buckling,  and  the  pre¬ 
buckling  displacements  are  assumed  to  be  small  so  that  the  additional 
displacements  can  be  related  to  the  undeformed  geometry  of  the  shell. 

The  additional  strains  are  taken  as  the  nonlinear  terms  of  equations  (11)  , 
namely*  , 


2 

x 


€’ 

9 


1  .2.2  2  _  .  „ 

— r - ~ —  (w,  +  v  cos  a  +  2vw,  cos  a) 

2x2  sin2  a  *  V 


(cor "inued) 


*These  nonlinear  terms  differ  from  those  used  in  previous  studies^* 10>25 
by  retention  of  the.  tangential  displacements  terms. 
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K  X  1U 
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(«S) 
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vies  tkt  mcI'jkjt  t 


cats  I  and  K  cf  cfutUu  C32J-(35). 

net  nsu  be  replaced  by  S  ♦  ■ 


xnl  j?  ♦  M  respectively,  • mi  the  potential  of 
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e  total  potential  energy  htriif  heck  1  leg  and  V  is 
the  strata  energy  stored  dvr  leg  hock  1  leg  fc»  the  adddle~serface  forces  of 
the  —berk  led  state;  Chat  is 


«.  *  v >;  ♦  •«  * 


(87) 


The  express loos  for 
.?«  chroegh  (M) . 


o  o 

ftr  axe  gives  bf 


Alter  epuiiou  (36s-(38»  with  the  aooliaear  tens  odttef 

iisj  surface  toad  terns  added  as  per  Appendix  I)  are  also  considered, 

tih-  stability  equations  and  boundary  conditints  are  obtained  as 
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L  * 

* 

2M 

X  ^  • 

V5?  \ 

•  *  ■ 

*  sir  xc.x 
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x  sin  a 


x  sin  a 


H  2M 

_2*5_  .  *9.0 

x  2 

x  sin  a 


cot  a  „  1  |  |, 

♦  — - - M  +  —  ’  |(«.  +  v  Co: 


-  r  sin  a  J 


:os  a)X 
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X  X  X 
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♦ «  I  I 


X  x® 
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x  sin  a 


K 


The  botaidar;  coni  it  ions  are* 

(a)  along  mt  and  x*x^ 


■  »  0 
x 


or  6u  -  0 


H  cus  a 

*  +  ;  ■ 
x®  x  sin  a 


or  6v  *  0 


sin  a  (xM)  -  M  sin  a  -  2K  (91 

XX  9  x®,® 

* 

+  x  sin  a  V,  I  +  (w,  +  v  cos  G0H _  -  0  or  6w  •  0 

x  xo  p  x»?o  — 


M  -  0 
x 


or  6w,  *  0 

—  x 


(b)  along  9  *  ®t  and  9  * 


H  -  0 
*9 


or  6u  *  0 


M  cos  a 


N  -  -*  -  «  0 

9  x  sin  Q 


or  6v  -  0 


M„.„ ' 2  sl“ a  <*%>  >«  <«> 
* 

+  (w,  +  v  cos  Q)N  +  x  sin  aw.  N  *0  or  dw  =  0 
9  90  x  *90  — 

(continued) 


Equations  (49)  -(51),  with  the  nonlinear  terms  omitted,  must  also 
be  used  in  the  derivation. 
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0 


or  6w.  *  0 

—  9 


M..  = 

v 

(c>  at  the  corners  of  the  segment 

M  =0  or  6w  -  0  (93) 

The  stability  equations  and  boundary  conditions  presented  above  may  also 
be  written  in  terms  of  either  forces  and  displacements  or  displacements 
only  by  using  equations  (39)  and  (40).  In  this  study,  however,  these 
rewritten  equations  may  be  obtained  directly  from  equations  (41) '(43), 
(46)-(48),  and  (52)-(54)  simply  by  deleting  the  inertia  terms  and  apply¬ 
ing  a  zero  subscript  to  the  forces  in  the  nonlinear  terms  consisting  of 
products  of  forces  and  displacements. 

The  stability  equations  and  boundary  conditions  derived  in  this  appendix 

represent  linear  buckling  theory  based  upon  the  strains  and  curvatures 

defineu  by  equations  (11)  and  (12)  respectively.  All  the  foregoing 

9  19  25 

equations  reduce  to  those  presented  by  Seide  ,  Singer  ,  and  Schnell 

* 

when  terms  marked  by  a  (  )  superscript  are  deleted. 


APPENDIX  III 


REDUCTION  OF  THE  GENERAL  NONLINEAR  CONICAL  SHELL  EQUATIONS  TO 
THOSE  FOR  A  CIRCULAR  CYLINDRICAL  SHELL 

In  the  present  study,  general  nonlinear  equations  for  the  bending, 

buckling,  and  vibrations  of  conical  shells  have  been  developed.  The 

geometry  of  a  truncated  conical  hell  if  such  that  it  develops  into  a 

circular  cylindrical  shell  of  radius  R  =  x  sin  O.  as  the  semi-vertex 

angle  a  approaches  zero.  Accordingly,  all  equations  developed  in  the 

present  study  for  a  conical  shell  reduce  to  those  for  a  circular  cylin- 

18 

drical  shell  (as  presented  by  Mayers  and  Wrenn  based  on  Sanders* 
theory)  through  the  simple  set  of  substitutions 

x  sin  a  =  R 


a  =  0 

1=0  (94) 

x 


where  y  is  the  tangential  coordinate  on  the  shell  middle  surface  in  the 
circumferential  direction.  Some  of  the  more  important  relationships 
derived  through  use  of  the  above  substitution  are  listed  in  this 
appendix  for  completeness. 

MIDSURFACE  STRAIN-DISPLACEMENT  RELATIONS  (rotation  about  the  normal 
neglected). 


By  virtue  of  equations  (94),  equations  (11)  reduce  to 
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2 

’x 


* 

•k 

€  =  v,  -  w/R  +  “(w,2  +  v2/R2  +  2vw,  /R) 

y  y  ^  y  y 
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7  =  v,  +  u,  +  w,  (w,  +  v/R) 
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(95) 
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CURVATURE-DISPLACEMENT  REIATICNS 


Equations  (12)  reduce  co 


K  -  v , 

X  XX 


H  *  W,  +  V,  /R 

cp  yy  y 

1  * 

x  *  w>  +  r  v,  /R 
xcp  xy  2  x 


Vfhen  the  terns  marked  by  a  superscript  (  )  are  deleted,  the  equations 
reduce  to  those  of  classical  von  Karman- Donne 11  theory. 


TOTAL  POTENTIAL  ENERGY 

The  total  potential  energy  is  given  by 

U,  *  U  +  U.  +  VT  +  V  + 

T  m  b  L  m  b 

Reduction  of  equations  (23)  and  (75)  through  (77),  in  conjunction 
with  equation  (97),  yields 

* 
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The  expression  given  by  equation  (98) ,  with  the  applied  load  terms 

18 

deleted,  was  used  by  Mayers  and  Wrenn  in  obtaining  a  first-approxima¬ 
tion  solution  to  the  free  vibration  behavior  of  a  circular  cylindrical 
shell  undergoing  large  displacements  and  moderate  rotations. 

FORCE -DISPLACEMENT  EQUATIONS  OF  MOTION 

Reduction  of  equations  (41) -(43),  with  the  addition  of  the 
applied  load  terms,  yields 

N  +  N  y  =  phii  -  q 

x, x  xy,  nx 

•jlf 

N  +  N  +  (D/R2)[R(w,*  +  w,  )  +  v,*  +i(l-v)v,  ] 

y, y  xy,x  L  xxy  yyy  yy  2  ’xxJ 


-  [(w,  +  v/R)N  +  w,  N  ]/R  =  phv  -  q 

l  v  >y  y  ’x  xy  y 


4  2 

Dv  w  +  D(V  v)  ,  /R  -  N  /R  -  [(w,  +  v/R)N  +  w,  N  ], 

y  y  y  xy  x  xJ  x 


-  t(w,y  +  v/R)Ny  +  w,xNxy],y  =  qz  -  phw 


BOUNDARY  CONDITIONS 


Equation  (49) -(51)  reduce  to 


(a)  along  x  =  x^  and  x  = 


N  =  N 
x  x 


or  6u  =  0 


N  +  M  /R  =  N 
xy  xy  xy 


or  6v  =  0 


M  -  2M  +  N  (w,  +  v/R) 

x,  x  xy,y  xy  y 


(100) 


+  w,  N  =  0  -  M 

xx  tc  xy,; 


or  5w  =  0 
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N 

y 


(105) 


The  compatibility  relationship  [equation  (61)]  reduces  to 


4_  2 

V  F  «  E(w, 

xy 


w,  w,  )  -  Ev,  / R 
’xx  ’yy  ’xx 


(106) 


and  the  lateral  equilibrium  equation  [equation  (62)]  becomes 


Dv  w  -  hF,  /R  -  h(w,  F,  -  2w,  F,  +  w,  F,  )  =  -phi#  (107) 
XX  '  ’xx  yy  ’xy  xy  ’yy  xx  K  '  ' 


Equations  (106)  and  (107)  represent  accurately  only  the  free  nonlinear 

vibrational  behavior  of  a  circular  cylindrical  shell  having  a  large 

number  of  circumferential  waves.  Various  approximate  solutions  for 

periodic  behavior  have  recently  been  obtained  for  these  equations, 

19  20  22 

notably  those  by  Chu  ,  Nowinski  ,  and  Evensen  .  The  most  recent 

solution,  based  on  the  energy  formulation  of  the  problem,  is  that 

18 

obtained  by  Mayers  and  Wrenn  .  In  their  work,  previous  periodic 
solutions  are  discussed  and  their  validity  is  questioned  on  a  minimum- 
energy-criterion  basis. 


Equations  (106)  and  (10”)  with  the  radial  inertia  term  omitted  are  one 

/  / 

of  the  equation  sets  commonly  referred  to  as  the  von  Karman-Donnell 
equations  governing  the  large  deflection  behavior  of  thin  plates  and 
shells  in  the  presence  of  midsurface  forces. 


All  of  the  conical  shell  equations  derived  in  Appendix  I  for  nonlinear 
behavior  under  arbitrary  loads  and  in  Appendix  II  for  stability  behavior 
may  also  be  easily  reduced  to  corresponding  equations  for  circular 
cylindrical  shells  through  the  substitutions  given  by  equations  (94). 
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APPENDIX  IV 


REDUCTION  OF  THE  GENERAL  NONLINEAR  CONICAL  SHELL  EQUATIONS  TO  THOSE 

FOR  A  CIRCULAR  FIAT  PLATE 


Ln  Appendix  III  the  limiting  case  of  a  circular  cylindrical  shell  is 
considered  and  the  basic  equations  describing  its  nonlinear  behavior  are 
easily  derived  from  the  conical  shell  equations.  Another  limiting  case 
occurs  when  the  semi-vertex  angle  a  approaches  the  value  n/2.  In  this 
instance,  the  general  nonlinear  conical  shell  equations  reduce  to  those 
for  a  circular  flat  plate,  with  a  concentric  hole,  by  virtue  of  another 
simple  substitution;  that  is, 

a  =  «/2 

x  sin  a  =  x  -  r  (108) 

INPLANE  STRAIN-DISPLACEMENT  RELATIONS  (rotations  about  the  normal 

~  neglected) 


By  virtue  of  equations  (108),  equations  (11)  reduce  to 


e  =  (u  +  v,  )  /r  4-  w,  2 J22 

9  9  9 


<%  •  v)/r  + 


(109) 


CURVATURE-DISPLACEMENT  RELATIONS 
Equations  (12)  reduce  to 


K 

r 


w. 


rr 


K 

9 


+  w,r/r 


K  =  W,  /r  -  w,  /r 
rep  rep  ep 


(110) 
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TOTAL  POTENTIAL  ENERGY 


The  total  potential  energy  is  given  by 


U„  =  U  +  U,  4  V,  +  V  +  VL 
i  :n  b  L  in  b 


(111) 


Reduction  of  equations  (23)  and  (75) -(77)  in  conjunction  with 
equation  (111)  yields 

"t  '  7 7X  ff  |<U'r  +  5  "4>2  4  f<“  4  ''■?>/r  4  W4/2r2]2 

*  9  r 

+  2v(u  +  -•  w  2)[<u  +  v,  )/r  +  w,2/2r2] 

I  l  r  (£  c 0 

L  “  v  2.  ^ 

+  —  [v,r  +  (u,^  -  v)/r  +  w,rw>cp/r]  j  rdrdcp 

+  f  +  [W’cpcp/r2  +  W’r/I]2  +  2v  "’r^V^  +  W’r/r] 

cp  r 

+  2  (1-v)  (w,^/r  -  w,^/r2) “  |  rdrdcp  -  ff  (q^_u  +  q^v  +  q^w)  rdrdcp 

9  r 

-  f  [r(N  u  +  N  v)  ]2  dcp  -  f  [N  v  +  N  u]2  dr 

J  r  r9  r  J  9  r9  i 

cp  1  r  m. 

•/  tr<-V>r  +  %Vr  +  V>)2  d<P  '/  t-VV4 


cp  1  r 

-  92 

+  M  w,  /r  +  Q  w]  dr 

r9  9  9  q, 

FORCE-DISPLACEMENT  EQUATIONS  OF  MOTION 

Reduction  of  equations  (41)  -(43),  with  the  addition  of  the 
applied  load  terms,  yields 

N  +(N  -N)/r  +  N  /r  =  phii  -  q 

r,r  r  cp  rcp,cp  r 


(112) 


(continued) 
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(N  +  2N  )  /r  +  N  -  phv  -  q 
cp.e?  r'/  r(f,r  > 


(113) 


where 


Dv‘  w  -  (l/r)[(w,  N  +  rw,  N  ),  +  (w,  N  /r  +  w,  N  ) ,  ] 

r  cp  rep  rrr  ep  cp  r  rep  cp 

=  -phi)  +  q 


'r  <  >  '  [''r  <  >i  -  4  t  (  >.tr  +  <  >.,/*  +  (  ),  /r2J 


BOUNDARY  CONDITIONS 


Equations  (49)  -(51)  reduce  to 

(a)  along  r  *  r.  and  r  =  r^ 


N  -  N 
r  r 


N  =  N 
rep  rep 

M  +  rM  -  M  -  2M 
r  r,r  ep  rcp,cp 

+  N  w,  +  rw,  N  *  rQ  -M 
rep  cp  r  r  r  rcp,cp 


M  *  M 
r  r 


(b)  along  cp  =  cp^  and  ep  =  cp2 


N  =  N 
rep  rep 

N  =  N 
cp  ep 

M  /r  -  2M  /r  -  2M 
cp,cp  T9  rcp,r 


(114) 


or  6u  *  0 
or  6v  =  0 


(115) 


or  6w  =  0 
or  6w,r  =  0 

or  6u  =  0 
or  6v  =  0 


(116) 


+w,  N/r+w,  N=Q-M  . 

cp  cp  r  rep  cp  rep,r  or  ow  =  0 


M  =  M 

cp  cp 


or  6w,  =  0 

—  cp 
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Equations  (58)  reduce  to 


rN  +H  •  H  +«  *  0 

r,r  r  cp  rep.cp 

N  +  2N  +  rN  =  0 
cp»cc  rep  rcp.r 


(118) 


DV  w  -  (1/r) [w,  N  +  rw,  N  )  ,  +  (w,  N  /r  +  w,  N  )  ,  ]  =  -phw 

r  rep  r  r  r  cp  cp  r  rep  cp 

Using  the  first  two  of  equations  (11S),  the  third  may  also  be  written  as 


D7«-  w  “  Nm(w»x  +  W’wyr2)  ‘  w*rr“x  +  2liJw*  J*2  ‘  w’  /r)  =  “phi; 
*  cp  x  epep  rr  x  xcp  dp  ^rp 


The  stress  function  definition  of  equations  (59)  reduces  to 


(119) 


N  /h  -  F  /r  +  F  ,  /r 
r  r  epep 

=  F. 


(120) 


N  /h  =  F,  /r  -  F,  /r 
xcp  rep  cp 

The  compatibility  relationship  [equation  (61) ]  reduces  to 
4  2  2 

V  F  =  (E/r  ) [w,  -  w,  w,  -  (w,  /r) ,  -  rw,  w,  ] 

r  L  ’rep  rr  epep  cp  r  ’r  rrJ 


(121) 


and  the  lateral  equilibrium  equation  [equation  (61)]  reduces  to 

DvJ  w  -  (h/r2)[w,  F,  -  2w,  F,  +  w,  F,  +  2(w,  F,  ) ,/r 
r  rr  epep  rep  rep  epep  rr  cp  cp  r 


+  r  (w ,  F ,  ),  Jr  +  r  (w ,  F ,  )  ,  -  2w,  F,  /r  ]  «  -phw 
repr  rrr  cp  cp 


(122) 


Equations  (121)  and  (122)  describe  the  nonaxisymetric,  nonlinear 
vibrational  behavior  of  a  thin  circular  plate.  When  axial  symmetry 


APPENDIX  V 


SUGGESTED  NUMERICAL  SOLUTION  PROCEDURE  WITH  ILLUSTRATIVE 
APPLICATION  TO  THE  NONLINEAR  VIBRATIONS  OF  A  BEAM 

Additional  definition  of  symbols: 


A 

I 

l,m*i,r 

L 


u 

w 

X 

z 


p 

$t(x) 


cross-sectional  area  of  beam 

cross-sectional  area  moment  of  inertia 

integers 

length  of  beam 

generalized  mass  functions 

generalized  coordinates 

axial  displacement  of  beam 

transverse  displacement  of  beam 

axial  coordinate  of  beam 

transverse  coordinate  of  beam 

mass  density  of  beam 

functions  of  x 

functions  defined  by  equations  (130) 


The  proposed  numerical  solution  for  the  nonlinear  vibrational  behavior 
of  a  conical  shell  can  best  be  explained  by  demonstrating  its  use  for 
the  simple  case  of  the  nonlinear  vibration  of  a  beam. 

The  expression  for  the  total  strain  energy  of  a  beam  in  the  absence  of 
applied  loads  may  be  obtained  directly  from  its  conical  shell  counter¬ 
part  [equation  (23)]  as 

L  2  1* 

u  =  \  y*EA(u,x  +  |  w,x)2dx  +  \  dx  (124) 

o  o 

and  the  kinetic  energy  is 

L 

T  =  j  J pA(u 2  +  w2)  dx  (125) 

o 

General  displacement  functions  may  be  assumed  in  the  form 
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w 


u 


l 

*  E  *<(*> 

i-i  1 
m 

-  E  <MX> 

i*Afl 


q^t) 

q1(t) 


(126) 


Substitu.ion  of  equations  (126)  into  equations  (124)  and  (125)  and  the 
enforcement  of  the  simultaneous  vanishing  of  the  first  variation  of  the 
Lagrangian  (T-U)  with  respect  to  the  various  q^(t)  lead  to  a  set 
of  m  second-order  nonlinear  differential  equations;  that  Is, 


E  M  +  *  =  0 
til  ri  1  r 


m 


E  Mrl  +  *r 
1=/+!  r 


1  SrSi 


jfrfl  s  r  s  m 


(127) 


where 


U 

Mri  *y*pA  ^(x)  A.(x)  dx 


(128) 


and 


!r"  /“[  £  *:WiHf  E  * 

r  o  1  i-Jfl  1  Jl  1=1 


$  ^(x)  dx 


Lt  ^ 

+  \  J  ea\  £  *  i<x)qt(t) 

O  Is  1 

+  J  ei*;  E  #  ,^(x)qi(t)dx 


$  ' (x)  dx 
r 


(Ursi) 


f  r  ^ 

r  -  \  J  EA  £  $  '(xjq^t) 

o  ^  1=1 


(129) 


<&  ^(x)  dx 


(continued) 
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+  f  EA$  '  (x)  $  '  (x)q  (t)dx  (i+1  s.  r  *  a) 

^  r  i-i+1  1  1 

o 

The  two  nonlinear  second-order  differential  equations  given  by  equations 
(127)  may  be  reduced  to  2m  first-order  nonlinear  equations  by  the 
substitutions 

^  =  Qi  i  =  1,  2 . m  (130) 

Equations  (127)  then  become 


A 

E  Mri  Qt  +  K  “  0 
i«l  ri  i  r 


E  M  Q  +  *  -  o 

i-j&fl ri  1  r 


1  £  r  £  l 


t+1  s  r  H  m 


(131) 


Before  the  numerical  solution  is  obtained,  functions  ^(x),  i*l,  2,  ... 
..m  must  be  assumed.  These  may  be  trigonometric  or  polynomial  in  nature, 
or  in  fact,  any  type  of  function  which  satisfies  the  geometric  boundary 
conditions  of  the  problem.  Once  these  functions  are  selected,  the 
problem  is  given  initial  values 

’i<0)  ■  v 

(132) 

q  (°)  - 

1 


and  the  solutions  are  obtained  numerically  using  a  greatly  modified 

30 

Runga-Kutta  technique  introduced  by  Rosser  and  slightly  modified 
■>1 

further  by  Skappel"'  . 

a 

For  a  given  time  increment  Vt,  the  initial  values  of  q^  and  q^ 
are  used  and  the  ijr  functions  are  computed  through  numerical 
integration  over  x  per  iteration.  This  process  continues  until  the 
required  history  is  obtained.  The  order  of  approximation  l  and  m  must, 
of  course,  be  given. 
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The  main  merit  of  this  proposed  numerical  approach  is  that  there 

is  no  need  to  expand  the  assumed  displacement  functions  by 

multiplying  out  the  infinite  series  and  collecting  terms  in  the  various 

's.  This  problem  is  prohibitively  severe,  especially  in  the  case  of 

nonlinear  problems  wherein  infinite  series  must  be  raised  to  the  third 

power  such  as  evidenced  in  equation  (129),  This  is  precisely 

why  the  bulk  of  nonlinear  structural  analyses  to  date  are  limited  to 

low-order  approximations;  in  the  case  of  the  nonlinear  vibrations  of 

thin  shells,  such  low-order  approximations  do  not  allow  for  the  devel- 

18 

opment  of  nonperiodic  motion  as  demonstrated  by  Hayers  and  Wrenn  .  It 
is  also  evident  from  equations  (129)  that  the  method  is  not  limited  to 
homogeneous  beams,  since  the  numerical  process  need  not  discern  any 
type  of  continuity  in  the  axial  direction. 

With  the  use  of  the  numerical  approach  proposed  here,  the  order  of  the 
approximation  (that  is,  the  manber  of  terns  retained  in  the  series)  may 
be  selected  at  will  and  limited  convergence  tests  may  be  run,  the 
limitation  being  the  amount  and  expense  of  computer  time  budgeted. 
Finally,  the  nuserlcal  process  itself  is  known  to  be  smooth,  and  no 
numerical  roundoff  difficulties  are  anticipated. 
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